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Abstract 

The b quark low-scale running mass w-kin is determined from an analy- 
sis of the T sum rules in the next-to-next-to- leading order (NNLO). It is 
demonstrated that using this mass one can significantly improve the conver- 
gence of the perturbation series for the spectral density moments. We obtain 
"^kin(lGeV) = 4.56 ± 0.06 GeV. Using this result we derive the value of the 
MS mass in: m.(fh) = 4.20 it 0.1 GeV. Contrary to the low-scale running 
mass, the pole mass of the b quark cannot be reliably determined from the 
sum rules. As a byproduct of our study we find the NNLO analytical expres- 
sion for the cross section e+e" — > QQ of the quark antiquark pair production 
in the threshold region, as well as the energy levels and the wave functions at 
the origin for the ^6*3 bound states of QQ. 

PACS: 11.65.Fy; 12.38. Bx,12.38.Cy, 11. 55Hx 



I. INTRODUCTION 



The value of the bottom quark mass is an essential ingredient of the theoretical descrip- 
tion of 6-hadrons. Among various applications, probably the most important one at present 
is the determination of the Cabibbo-Kobayashi-Maskawa matrix elements from 5-decays. 
Determination of the b quark mass is based on the sum rules for the T mesons, proposed 
about 20 years ago [jl],|^ . In the past several years, the sum rule analysis has been undertaken 
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several times by different authors. In particular, in the paper a very high accuracy of the 
b quark pole mass was quoted. The next attempt to extract the precise value of the b quark 
pole mass from the sum rules was undertaken by Jamin and Pich in [Q. Their result differed 
from that of Ref . . The origin of the discrepancy between these two results, as well as the 
flaws in both derivations, were pointed out by Kiihn et al. 0, who also determined the b 
quark mass. Two most recent papers on the subject were devoted to an improvement 
of the theoretical accuracy of the mass determination and to a more realistic estimate of the 
theoretical error. 

In parallel to these developments, it became more and more clear in the past years that 
the concept of the pole mass of a heavy quark is not a good one due to the intrinsic ambiguity 
of the order of Aqcd in its numerical value P,[lO|. In contrast to this observation, all previous 
analyses were aimed at determining the pole mass of the b quark from the sum rules. We 
note in this respect that practical problems in attempts to use the pole mass in heavy quark 
physics are well appreciated; one of the vivid examples is provided by the calculation of the 
inclusive semileptonic decay widths of the B mesons (see for a review). In response to 



this problem, it was pointed out how the "proper" quark mass relevant for non-relativistic 
problems can be defined |T^. In this paper we try to determine this properly defined mass 



from the sum rules. 

For the theoretical analysis the moments of the photon polarization operator are used. 
These moments can be computed analytically and compared to the experimental ones. The 
moments of the photon vacuum polarization function are defined through the dispersion 
integral: 



where R{s) is 



a(e+e ^bb) 47ra^ED("^b) 
K{s) = , CTp = — . [2) 
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We defined the moments to be dimensionless by multiplying them by the mass Mi of the 
first T resonance in a suitable power. Also, = — 1/3 is the electric charge of the b quark 
in units of the positron charge. 

The moments Ain can be calculated using experimental input for R{s). One gets: 



Mf" / Qtt ^ Tfc f Rc{s)ds 



oo 



where and are the masses and the electronic decay widths of the first six T resonances 
and aQEol"^?') — ^■^'''^qed used. The Rc{s) describes the experimental spectral density, 
associated with the energy region above the open BB threshold; it is rather poorly known. 
The necessity to suppress the contribution of this region is one of the reasons for using as 
high values of n as possible. 

On the other hand, the same moments can be calculated theoretically, using the Operator 
Product Expansion (OPE) for the photon polarization function. The first non-perturbative 
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correction to the moments is associated with the ghion condensate [0,0. It was shown that 
this contribution grows with n; however, for n < 20, the non-perturbative contribution to the 
moments was estimated to be less than one per cent 0,^. Given exponential sensitivity of the 
moments to the value of the h quark mass, the influence of the non-perturbative corrections 
is clearly minor and can be neglected. Therefore, the whole analysis for n ~ 10 reduces to a 
careful treatment of the perturbative effects in R{s). However, the perturbative treatment 
is not simple, since for large values of n the dominant contribution to the perturbative 
moments comes from the threshold energy region. The relative velocity of the hb system 
there is of the order of a^. In this case, the theoretical spectral density can be calculated 
in the framework of the non-relativistic QCD, which means a simultaneous expansion of 
the spectral density in as and in the relative velocity (3 of the quark antiquark pair. As 
is well known, the standard perturbation theory is not adequate in the threshold region 
and the leading order approximation is the solution of the Coulomb problem, which resums 
all corrections of the form [ag/ (3)^. Going to NLO and NNLO, one calculates the spectral 
density R{s) in the threshold region resumming all 0{{as/ P)^ x [1; a^, /?; a^, asP, P^]) terms. 
For this purpose we use the so called direct matching procedure which is described e.g. 
in p^ , p!5| ,|7| . We will not discuss all necessary details of this approach here and will merely 



quote the results of the calculations. On the other hand, a part of the NNLO corrections 
was treated numerically in |]l^|15l; for this reason, we present some additional theoretical 



results, which provide the imaginary part of the polarization operator in the threshold region 
in completely analytical form to NNLO. 

The theoretical expression for the spectral density employed in this paper reads: 



R{s) = limlm 







(l-|^)G(r,0) 




(4) 



where G{r, 0) is the Green function of the non-relativistic Schrodinger equation 

{H-E-i6)G{r,ri)=6^^\r~ri), E = ^s-2m. (5) 

Taken literally, Eq. (^) is ill-defined due to improper treatment of the limit r — > within 
the non-relativistic approach. In what follows, we circumvent this difficulty incorporating 
the full QCD result and matching it with its non-relativistic counterpart. 

Integrating the theoretical expression for the spectral density in Eq.(|I|) and equating the 
obtained result to the experimental moments Eq.(|^) one obtains the sum rules of the form: 

= M'^'"', (6) 

which give the numerical value for the b quark mass, potentially with a small uncertainty. 

To summarize, the usual (and commonly used for the analysis) statement about theo- 
retical calculations of the large-n moments can be formulated as follows: 1) for n ~ 10 — 20 
the OPE ensures that non-perturbative corrections are small; 2) for such n the perturba- 
tive calculation of the moments requires the resummation of the Coulomb enhanced terms 
0{asy/n). This is achieved by calculating the spectral density R{s) in the threshold re- 
gion which requires a resummation of the x [l;as,P;al,asP,P^]) terms for the 
NNLO accuracy. It is often assumed implicitly, that this picture is correct independently 



3 



of all other parameters entering the sum rules analysis. We do not believe that this is the 
case. In particular, for the first point to be correct, it is crucial that one does not attempt 
to extract the pole mass of the heavy quark from the sum rules. 

When this paper was prepared for publication, two other papers on this subject appeared 
I^J^ where the NNLO analysis of the T sum rules has been performed. The aim of both 
papers was to extract the pole mass of the b quark. In our approach to the same problem we 
treat the so called low energy running mass as a quantity which can be determined within 
the sum rule analysis; the pole mass of the b quark is used at the intermediate stages of the 
calculation only. 

Let us present some heuristic arguments in favor of this approach. It is known that 
the pole mass of the quark cannot be defined when non-perturbative effects are addressed 
PJTO[]. It is also believed that the bad behavior of the perturbation series in the relation, 
say between the pole mass and the MS mass, signals this. The consequence of these facts 
is that there is an irreducible ambiguity of the order of Aqcd in the numerical value of the 
pole mass. If then the pole mass of the quark is used in the sum rules analysis, its infrared 
sensitivity leads to new infrared effects, which have no counterpart in the standard OPE . 
In particular, they are not described by the gluon condensate [pH] . 

On the other hand, one can realize, that this is an artifact of the adopted procedure and 
the easiest way out of this problem is to abandon determination of the pole mass from the 
sum rule analysis. Therefore, we use the pole mass only as a tool to write the expression for 
the non-relativistic Hamiltonian; however, we do not treat the pole mass as a fixed number 
and recalculate it consistently, order by order in perturbation theory. 

Instead of the pole mass, one should determine some "proper" mass, which does not 
suffer from a numerical ambiguity due to contributions of the soft momenta region. Such 
proper masses are known - one of the most familiar in this respect is the MS mass. This 
will not be our choice, however. 

Elementary physical considerations suggest that the threshold problems are the low-scale 
problems, in principle. The typical scale is /x ~ m/ \fn ~ 1 — 2 GeV for ~ 10. A useful and 
reliable mass should therefore be normalized at such low scale -C m. On the other hand, 
in order to suppress the contribution of the infra-red region, the inequahty Aqcd ^ should 
be respected. The use of such low-scale running masses in various aspects of heavy quark 
physics was repeatedly advocated in the last years (for a detailed discussion and further 
references see [|lT]). 



The low-scale running mass cannot be defined uniquely, because the only purpose of 
such definition is to remove the uncontrollable contribution of the soft momenta region 
which affects the pole mass. In this paper, we will work with the so called kinetic mass 
suggested in W, 



The relation between the pole and the kinetic masses is known to second order in |16 



mpole = mkin(/iQ) + [A(/iQ)]p,rt + 2mkin(/iQ) V'-^^'^>\ pert ' 

where 



3 ' TT 1 TT 



= IC.yH^ \ 1 + 
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pert 



2"s(/ily 



1 + 



TT 



TC 



12 2 yuJ^ I 6 12 



For the rest of this paper, denotes the strong couphng constant in the MS scheme. 
Also, we choose /iq to be equal to 1 GeV, which seems to be a reasonable choice for the 
problem at hand. Then the ratio /ig/m is of the order of and this gives the counting rule 
for the contributions to the mass which should be accounted for when one goes from one 
order of perturbation theory to the other. For example, to obtain the LO result we consider 
the pure Coulomb potential without any corrections in the non-relativistic Hamiltonian for 
heavy quark antiquark pair. Correspondingly, the LO relation between the pole and the 
kinetic masses is: 



'^pole 



a. 



?^km(/iQ) + -Cp—^lQ. 
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The NLO and the NNLO corrections to this expression are added in accordance with the 
above counting rule. Also, for our treatment, we use the same normalization scale /ii for the 
strong coupling constant in the expression for the mass, as is used in the non-relativistic 
Hamiltonian. This scale is called /isoft in the rest of the paper. 

We would like to note that, according to the above counting rules, it would be necessary 
to know the term a^fj^q in the relation between the pole mass and the kinetic mass. This 
term is not known at present. The BLM-type estimate of this term is however available 
and can be extracted from [M. Explicitly, the necessary term reads: 



6 [A(/.)] 



pert 



In' 




(8) 



where /Sq and (3i are given in the next Section and c/i 2 read: 




2 r 




In 2 - 



(9) 



In the numerical analysis of the last Section, we check the sensitivity of our results to the 
possible modification of the d2 term due to additional terms which are not accounted for 
in the BLM approximation. We find, that our final result for the kinetic mass is rather 
insensitive to it. 

As the result of our analysis, we find that the perturbation theory for the pole mass is 
not applicable: typically, the NNLO corrections to the pole mass exceed the NLO ones and 
the dependence of the result on the choice of the scale of the strong coupling constant is 
very strong. These unwelcome features, therefore, do not permit a reliable determination 
of the b quark pole mass from the sum rules, with a trustworthy estimate of the theoretical 
uncertainty. 

On the contrary, the situation with the low-scale running mass looks more healthy: the 
perturbation series seem to be sign alternating and the dependence on the normalization 
scale for the coupling constant appears to be reduced, as compared to the pole mass. 
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The rest of the paper is organized as follows: in the next Section we discuss the framework 
of the calculation. In Section 3 the corrections to the Green function due to corrections to 
the static quark antiquark potential are derived. In Section 4 the corrections to the Green 
function due to relativistic corrections to the heavy quark Hamiltonian are obtained. In 
Section 5 we combine these results and present the NNLO expression for the theoretical 
spectral density R{s) in the threshold region. In Section 6 the results for the energy levels 
and the wave functions at the origin for the ^S^ QQ resonances are derived. In Section 7 we 
present our final analysis for the sum rules and determine the low-scale mass of the b quark. 
Finally we present our conclusions. 



II. THE FRAMEWORK OF THE CALCULATION 



We first discuss a framework of our calculations and introduce all relevant notations. As 
we mentioned already, in order to obtain the expression for the theoretical spectral density, 
we have to calculate the expression for the imaginary part of the polarization operator in 
the threshold region. The threshold region is characterized by a small value of the quark 
velocity /3: 



(3 



Aw? 



< 1. 



(10) 



Here and below m is the pole mass. The pole mass enters the usual perturbative expansion 
in quantum mechanics for the nonrelativistic quarks. Later we will extract the low-scale 
running mass from the pole one. 

Dynamics of slow moving quark antiquark pair is governed by the non-relativistic Hamil- 
tonian: 

H = H^ + V^{r) + U{p,v), 
riQ — , 



VAr) 



m r 



2/?oln(2/xV) + ai 



^7(p,r) = 



+ ( [Po (^41n2(/.V) + - ) + 2(A + 2/3oai) ln(/.V) + a, 

^ TrCpas 3, , Cpcis 
+ —^i^) 



+ 



4m^ 

3C f(^s 



SL- 



CpCls 



2m^7 



3 ^ ' ^ ' 2mr2 



In the above equations, the strong coupling constant is evaluated at the scale //soft- 

as = Q!s(/^soft)- 

The scale equals to fie^'^. 



(11) 



(12) 
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The operator U{p, r) is the QCD generahzation of the standard Breit potential [jlH . 



The last term in the expression for the operator U{p, r) is the non-Abelian contribution, 



originating from a correction to the Coulomb exchange, caused by a transverse gluon [jT9 
The potential Vi(r) represents a deviation of the static QCD potential from the Coulomb 
one. It was calculated to order in and to order in The coefficients there read 
explicitly: 

M 2(1 

Pi = -C\ - -CaTrNl - 4CfTrNl, 

ai = —Ca- —TrNl, 
y y 



/4343 ^2 ^ 22 . ^2 
V 162 4 3 I ^ 

1798 56 /55 /20. 



+ yCa) C^T^iVL - (y - leCa) C^T^iV^ + (y^^^^) ■ (13) 

The SU{3) color factors are Ca = S,Cf = 4/3, Tr = 1/2. A^'^, = 4 is the number of 
quarks whose masses have been neglected. 

Given the Hamiltonian H, one can find the Green function of the Schrodinger equation: 

{H-E-i6)G{r,ri) = 6^^\r~ri), E = 4l-2m. (14) 

Once the Green function is found, the non-relativistic cross section for the QQ pair produc- 
tion in e^e~ annihilation is obtained using Eq.(^. 

Treating the corrections to the Green function in the perturbation theory, one can con- 
sider the corrections due to Vi(r) and ?7(p, r) separately. The corrections to the Green 
function due to [/(p, r) were recently calculated in [|14|,|T5|. These corrections are not sim- 
ple conceptually, because they deliver divergent contributions to the Green function at the 
origin. The divergences are removed by matching the result of the calculations in quantum 
mechanics to the result of the full QCD calculation Technically, however, the calcula- 
tion of the correction caused by t/(p, r) is very simple and can be performed algebraically 
(see |]T3[). On the other hand, the corrections to the imaginary part of the Green function 
due to the V\ ir) perturbation can be calculated within the quantum mechanics and for this 
reason these corrections are rather simple conceptually. However, they provide the most 
challenging part of the whole calculation from the technical viewpoint. For this reason, the 
calculation of the corrections to the imaginary part of the Green function due to the V\{r^ 
perturbation is discussed below in some detail. 

III. CORRECTIONS TO THE GREEN FUNCTION DUE TO Vx{R) 

A. The Coulomb Green function 

In this section we collect useful formulas for the Coulomb Green function. A convenient 
expression for the S'-wave Coulomb Green function can be found in [|23|: 
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G(r,n) ^ J_|L= (i±l)%-(-..V.(4.v??Tv^). (15) 

In the above expression, u = CFas/{2j3) and k = m(3. 

Using this representation, one easily obtains the expression for G(r, 0): 



-imk ji^j. f f ^ + tY'^ 2ikrt 



G{r, 0) = — e*'^" / dt — e''"''. (16) 



From the expression for the cross section (cf. Eq. (^)), it is clear, that one is interested 
in the behavior of the Coulomb Green function for small r. For small r, the Green function 
diverges like 1/r; the principal divergence is related to the behavior of the free (a^ = 0) 
Green function: 

G^°nr,0) = £:e- (17) 



B. Generating function 

To calculate corrections to the Coulomb Green function at the origin caused by the Vi(r) 
perturbation, it is convenient to introduce a generating function g{cr): 

g{a) = JdhG\r,0)^^^. (18) 

Once this function is found, one easily obtains a correction to the Green function at the origin 
caused by the Vi(r) term in the potential. For further calculations, it will be convenient to 
separate the generating function into two terms: 

fl'(0") = Qfreeicr) + 5'l('^), (19) 

where 

g,{a) = J dh {G\r, 0) - G^(r, 0)} (20) 

In the above equations, Go{r, 0) is the free Green function (cf. Eq. ([T7|)). 
It is an easy task to calculate Qfree- One gets: 



We need further gi{(j). For the Coulomb Green function we use the representation from 
Eq. ([161) . Integrating then over r in Eq. (pO]), one gets: 
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An \ k 



CO 

r(2 + cr) J dtds 



l + tV" fl + s 



il + t + s) 



-2-0- 



(22) 



To proceed further, it is convenient to introduce new integration variables, 



t 



T 



Then one gets: 



1 + 



1 



l + s 



(23) 



2-0- 



Substituting also rj = pr we arrive at: 

,,.H^(|yrp..)/.,^/=^a-.r(:-^ 

^ r/ 

Finally, changing the variables t — > ^ with ^ = (1 — r)/(l — 77), we get 
^'^^'^ r(2 + a) / rfr/ J^Ji^. / ^^^^(1 - (1 - (1 - 



(24) 



(25) 



4:71 \k 



(26) 



The last integral is proportional to F2i(l + a, 1 + a; 2 + 2(t; 1 — ?]) and equals to 



(1- (1-^)0 



n=0 



(27) 



where the series representation for the hypergeometric function was used. Here {z) 
r{z + n)/r{z) is the Pochhammer symbol. Integrating over rj, we obtain 

g,{a) = — (-^ j r(2 + ana) [2V^(n + 1) - 2V^(n + 1 + a) - 5„] x 



r(rz + 1 — iz/) 



r(n + i: 



[T{n + l~iu + a) r(ri + 1 + cr) 
which can be presented in a more compact formQ: 

m^/iayi + a^^ (T^in + a) r(n - iu) r(n + a)] 



(28) 



Air \ k J a "\ V^in) V{n — iv^a) r(^) 



(29) 



In Eqs.(28j29|) dn stands for the partial derivative with respect to n. Eqs. (T9|,21,29) provide 
a necessary expression for the generating function. 



^We thank O. Yakovlev for pointing this to us. 
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C. Correction to the Green function due to ln"(2/ir)/r, n < 2 perturbation 



To evaluate these corrections we will use the representation of the generating function 
g{a) given in the previous section. For the case of interest, the correction to the imaginary 
part of the Green function at the origin due to ln"(2//r)/r perturbation is obtained as the 
n-th derivative of the imaginary part of the generating function with respect to cr at cr — > 0. 
To obtain these derivatives, we first expand the generating function up to the second order 
in cr: 



9iiee(cr) = 



m 



1 



An 1 a 



In 



in 



lE + 



In 



In^ 



TT 



111 



lE 



TT 



7b) +-(ln^-7£) -2C; 



+ 



and 



where 



g'M 



^i(^) = ^(^i(o) + ^^i(o) + yyi'(o) + ..., 



n=l 

(in^ + l) ^,(o) + ^ + la>V'(i-i^) 

-j^ oo 
^ n=l 



(30) 



(31) 



(32) 



(33) 



^i(O) 



1 °° r 3 ' 

~q H 1 9 (SV^^H - Ail){n)il){n - iu) + ^^(n - iz/)) 

n=l ^ 

n — iv) — ij\n-iv)}. (34) 

Using these expressions, we easily find the corrections to the imaginary part of the 
Green function at the origin, caused by ln"(2/xr)/r perturbations. In the formulas below we 
disregard the 1/a pole which is present in the expression for g^free since it does not contribute 
to the imaginary part. We obtain explicitly: 



SGM^-JdhG{r,0) 



,ln(2/ir) 



-m 



In- -7. 



+ 



TT 



6 



+^;(o) 



r 47r ^ 

In a similar manner one obtains the correction due to ln^(2//r)/r perturbation: 

,2W(2//r) _ 



(35) 



SGl2{I^) 



d'^rG'(r,0)^ 



-m 



In --7. 



+ 



TT 



ln|-7i.)-2C3 



+ <(0) 



(36) 
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D. The second iteration of the logarithmic perturbation 

The leading term of the static potential Vi(r) provides an 0{as) correction. Therefore, 
one should calculate the second order correction induced by the perturbation ln(2/ir)/r. 
Such calculation is described below. 

Explicitly, we have to calculate: 

5G(2)(^)= [ dhd^riG{r,0)V{r)G{r,ri)V{ri)G{ri,0), 1/(r) = ln(2^r)/r. (37) 



For this purpose we use the representation for the Green function G(r, ri), given in Eq. (|T5|). 
We define: 

^ oo 

-(T)'" r(-.)r(-..) /;^^'^-^''- '''' 





where 

n 



F(r,ri) = / d=^rd=^riG(r,0)^ G(r,n)- G(ri,0). (39) 

J r r-i 



Then 



rj2 

= Gi(a,(Ti)U=o,.,=o. (40) 



Using Eqs. ([T5|JT6D for the Green functions, and performing the integrations over r, ri 
and the variable t which enters the integral representation of the G{r, ri) function, we arrive 
at the following result: 



F{r,n) = -— -Fi(r,n), (41) 

87rK(l — lu) 



^ dudv fl + uV" fl+vV" ^ ( . . (A-l)(Ai-l) 



where 

A = l + u + r, Ai = H-'j; + ri, (43) 

and F21 (a, 6; c; z) is the Gauss hypergeometric function. 
The function Gi((J, cti) can be written as 

where 
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"^•'-'-''- r(-.)r(-..) /7^^''"'^-'' '''' 

For further calculation, it is convenient to define new variables 

u ^ Ui = u — t; V ^ Vi = V — Ti 
and integrate over r and ti. Performing another variable transformation, 

Ui Vi 

ui^x = — vi^y = — , 

1 + Ml I +Vi 

we finally get: 

L [1 — lu — a)L [I — lu — ai) \l — x/ \l — y J 

F2i{-iv, -a; I - iv - a; x)F2i{-iv, -ai, 1 - iu - ai, y)F2i{2, 1 - iv;2 - iv; xy). 

One notes that if the last hypergeometric function in the above equation is expanded in 
Taylor series in xy, integrations over x and y factorize. These series would, however, diverge 
for xy = 1. Nevertheless, upon integration over x and ?/, one gets a series which converges 
as We conclude therefore, that this operation is legitimate. We write W{a,ai) in a 

factorized form: 

W{a, a,) = ^^^-'"11]-'''^ . £ -^Tin, a)T{n, a,), (46) 

r(l — iv — cr)r(l — lu — (Ti) n — IP 

where 

1 

T(n,a)= [ dxx''-^-''' (^^) " F2i{-iv, -a-l - iv - a;x). (47) 



X 



According to Eq.(|40|), one needs an expansion of T(?t,, a) in a up to the first power. This 
can be easily done by expanding the hypergeometric function in Taylor series in Eq.(^), 
evaluating resulting integrals, extracting the limit cr — and then resumming the resulting 
series. We obtain: 



where 



T(n,a) = ^--(7Ti(n), (48) 
n — IV 



^ (^\ - I i){\-iv) i){n - iv){n ^ iv) iv 

— — H \ 7 r^; \ 7 —r^- i^yj 

n — IV n n(n — iv) n(n — ivV 



With this result, we obtain the series representation for the function W{(t, ai) with the 
required accuracy in a, cxi. Upon differentiation over a and ai at cr = (Xi = we get the 
correction to the Green function 5G'^'^^: 

onk ~^ n — IV \ n — IV 
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IV. CORRECTIONS TO THE GREEN FUNCTION DUE TO THE U{p, r) 

PERTURBATION 



We now briefly discuss corrections to the Green function at the origin caused by the 
operator f/(p, r) from Eg . (pH]) . This correction is obtained as 

6G^ = - J d^rG{r, 0) U{p, r) G(r, 0). (50) 

As long as we are interested in the QQ pairs, produced in the triplet ^-states, only the 
corresponding projection of the operator f/(p, r) should be considered. Substituting = 1 
and SL = in Eq.(|Tl|), it is easy to get that f/(p, r) can be presented in the following form: 

At this stage, it is advantageous to express this operator in terms of the zeroth order 
Hamiltonian i^o in order to apply the equation of motion for the Green function G{r, 0): 

iHo-E)G{r,0) = 5^^\v). (52) 

This is most easily done using the following commutation relations: 

. 47r^(3)(r) , 2L^ C^a, 
[Ho, ipr] = + ^ (53) 

H„1\^Ih^'-^^^,^, (54) 

where Pr = —i{dr + 1/r) is the radial momentum operator. One finds that the operator 
U (p, r) from Eq. ( |5lD can be written as: 

t/(p r) = - [ho -] + [Ho ipA - (2^^ + 3C'A)CFag 

' 4m 4m I ' r J 12m ' ^ 6mr^ 

Let us consider the first three terms of Eq. (^). Inserting them into Eq. (|50|) and using 
the equation of motion for the Green function (p^, we find: 

r2 



G(r,0) + 2^p,G(r,0)+ / d3rG(0, r) |^ + ^%^| G(r, 0). (56) 

J 4m Imr 



E ^ 3CpCis 



2m 2mr 

All terms in the above equation which cannot contribute to the imaginary part of the Green 
function have been omitted. 

Also, it is easy to recognize that the terms in Eq. (^), which still have to be integrated 
over r, can be easily obtained if one redefines the eigenvalue and the coupling constant of 
the lowest order Hamiltonian Hq: 
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ff„^«.P!_^Z-(:,|£); E^S^E^f^^l^p^. (67) 

m r \ 2m ) Am m 1 — 

The new Hamiltonian Ti is still of the Coulomb form and, therefore, the solution for the 
Coulomb Green function presented in the previous section can be used. 

Therefore, we conclude, that the only non-trivial calculation required here is the correc- 
tion to the Green function at the origin caused by the perturbation, which is explicitly 
given by the last term in Eq. (|55|) . We refer the reader to the paper [|l^ where the details 
of the calculation are discussed. 



V. COMPLETE RESULT AND MATCHING 

We now combine the results of the above calculations and write the final result in the 
form: 



R{s) = /C(/ihard,/ifact) (-Rl(s) + R2{s)) 

Rlis) = ^N,QlCFaslm(^H{CFas,/3) 



(58) 



l + CFaU^ + ^)H{CFa,,P) 



m^ 



/5o + 7^(A + 2/5oai) 
47r 



6Gr, ( ^ 



2^3 



2 



4(3^CFai 
(47r)2 



Ai6 



(59) 



In the above expression, i?2 is the contribution due to the Vi(r) perturbation. The scales 
/xi, /ia, fib there read explicitly: 



AH = Atsoft exp [7^;] , Ata = /ii exp 



2(3o + f^2(A + 2/3oai) 



/ife = yUi exp 



ai 
Wo 



(60) 



The function if(a, j3) was first obtained in a similar context in ||13[ and is given by: 



H{a,f3) 



(61) 



We have absorbed all energy-independent divergent contributions to the factor 
A^(/^hard5 Z^fact), which is determined by matching the above result for R{s) to the result 
of the perturbative calculations in full QCD p2[ in the region <^ /? -C 1, where both 
results are supposed to be valid. One gets: 



ir{,, ,, ( (^sijMmidf] ^ ^ /^as(Rard) 
/^lA^hard, /^factj — i + ^l^F ) + ^2^F 



TT 



TT 



(62) 



Ci = -4; C2 = CfC^ + CaC^^ + TrNlC^ + ThNhC^ - Ci ^ In f ^ V (63) 

4 VA'-hard/ 



14 



and 



39 



4 



^^^ln2-^ 
18 



TT^ , VP? 

— m —2 — : 

3 A'-fact 



151 13 



36 2 
44 4 



/179 8 
172" ~ 3 



In 2 



TT 

9 9 
11 



— m —2 — : 

2 A'-fact 



(64) 



VI. CORRECTIONS TO THE ENERGY LEVELS AND THE WAVE FUNCTIONS 

It is known, that the proper expression for the Green function which is vahd in the whole 
threshold energy region is given by the expression: 

G{E + %e- 0, 0) = y + / — . (65) 

^ ' ' ^ ^En-E-ie 27iEk-E-ie ^ ^ 

Here En and are the energy levels and the wave functions at the origin of the perturbative 
^5*1 QQ resonances, which can be calculated order by order in perturbation theory. On the 
other hand, when a correction to the Green function is calculated as a power series over a 
perturbation, the energy denominators entering the exact expression ( |65D are also expanded 
in powers of a^. It is therefore possible to extract the corrections to the energy levels 
and to the wave functions of the resonances by performing the Laurent expansion around 
E = E^^ of the corrections to the Green function obtained in the previous sections. On one 
hand, these results are interesting by itself, for applications to various problems that involve 
perturbative calculations for bound states both in QED and QCD. On the other hand, they 
are used below to construct the proper theoretical expression for the large-n moments. 
In the formulas below we denote: 

i^z = 'ipiz), = -^^{z), etc. 
To present our results, it turns out to be useful to define a function 

n-l 

k 



S.{n) = E 1^- (66) 

k=l 

The energy levels and the wave functions at the lowest order are given by: 
^(0) _ m{CFas{f^)Y hj,(o)|2 _ (mCjra,(/i))^ 

First we present an expression for the energy levels valid up to (relative) order (9(a^): 
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+ 



27r 



(L(/ii) + ^Z-^+i) - ^Z-^+i - 2 + — 



n 



2 r 



+ - 1) - 1 - - nV^„ + - {i)n+i + 7i?) 





2 




V16n 


~ 3 ~ 


1)} 



(68) 



where = /^soft, 



L(/i) = In 



/ /in 



and the scales are defined in Eq.(|60l). This resuh is in agreement with the one 



obtained in [El 



Using the same notations, we obtain the result for the square of the wave functions at 
the origin: 

|V^.|2= /C(/ihard,/ifact) (^) (/^O + ^ (A + 2/?o«l)) [3 (L(/^J + V'.+l) - 



2 \ 27r 



(T( \ ^ I ^2 ;/ +7S , 

n 3 



27r 



I (2nV^; + 27/^„ + 27^; + 1) + 2n< (t/-, - 7^;) + (^1 - 27^; - ^) 



+2nS2{n + 1) - 2nC3 + 2(2(1 + Ieu) + - - 27I 



n 



+ 



27r 



(4n2V;; - 2nV^„ + n(5 + 47^) - 6) L{^i,) 



n 



+ — ^ + ^ + i^nin + 2n(l + 37i=;) - 5) - 

ipninipn + 3(n - 1) + ^lEn) (6 - 5n + + 272^(3 - ^^(4) 



+ 



n 



27^(2^3(3 -n^Cs- I + n) 



n 
37 



+ 472^^3 (n) -2n52(^) 



+ 



24„^ vl+§)r(iS^)+'"»+'"=-i 



(69) 



In the above expression, /C(/ihard; /^fact) stands for the NNLO hard renormalization factor 
given explicitly in Eq.(|62|). 



16 



VII. THEORETICAL MOMENTS AND NUMERICAL ANALYSIS OF SUM 

RULES 



The theoretical moments can be conveniently separated into the contributions of the 
perturbative resonances and of the perturbative continuum: 

Mn =Vn+ Cn. (70) 

The resonance contribution reads: 

^^^^ k=i m? [l + 

and the continuum contribution is defined as the integral of the function R{s) / s^^^ over s 
(see Eq.(^)) taken above the threshold. 

We stress that the above expression for the moments differs from the result one gets, 
merely integrating the corrections to the Green function. The difference is due to the fact, 
that in Eg. ([TTD , we calculate the corrections separately to the numerator and the denomina- 
tor. Working in the limit ^/nas ~ 1, it is possible to expand the denominators in Eq. (|7ll) 
around their values for the exact Coulomb problem; this would produce (parametrically) 
an 0{as) corrections to the moments. However, there is a serious numerical difference be- 
tween the expanded and not expanded denominators. The origin of the problem is related 
to the large values of and ai^2 entering the potential. This effectively translates into the 
large values of the corrections to the energy levels; for this reason, the expansion of the 
denominators is not justified in our opinion. 

In contrast to the contribution of the perturbative resonances to the moments of the 
spectral function, the contribution of the continuous spectrum to the moments behaves 
nicely, as far as its perturbative expansion is concerned. The contribution of the perturbative 
continuum to the theoretical moments is obtained by numerical integration. 

For numerical analysis of the sum rules, we use the value of the strong coupling constant 
as[Mz) equal to the world average value as{Mz) = 0.118. This value of the strong coupling 
constant is evolved down to a required scale n using the two-loop renormalization group 
evolution equation. We will later comment on the sensitivity of our results to the value of 
the strong coupling constant at the Z-resonance. 

We also parameterize the unknown contribution of the experimental continuum Rc{s) in 
Eq.(^ by a constant, which we vary between and 2. The value of sq in Eq.(^) equals to 
(2 X 5.927 GeV)^, i.e. the continuum contribution starts at the threshold of the open BB 
production. To suppress the influence of this unknown contribution we have to go to rather 
high values of n. We have chosen n equal to 14, 16, 18, for our analysis. 

We fix the values of the hard and factorization scales at 5 GeV and examine the value 
of the 6-quark mass as a function of the soft scale /Xsoft- Our results are presented in Table 
1. The values of the kinetic mass at /xq = 1 GeV are found using Eq.(|^) treated with the 
necessary accuracy. 

The last column in Table 1 demonstrates how the perturbation theory works for the 
theoretical moment Aiu at LO, NLO and NNLO if the NNLO value of the kinetic mass is 
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Scales 


Order of PT 


14 


16 


18 


PT for M^A 
for NNLO mkin 




^kin 




^kin 






//soft =4.5 GeV 

A^hard = 5 GcV 

UfaM = 5 GcV 


LO 


4 fiQ 


4 57 


4 7 


4 57 


4 7 


4 575 


0.97 


NLO 


4.74 


4.49 


4.75 


4.50 


4.76 


4.51 


0.61 


NNLO 


4.89 


4.51 


4.90 


4.52 


4.91 


4.53 


0.69 


Hsoit = 3.5 GcV 

/^hard = 5 GeV 

Ufo^i ~ 5 GeV 


LO 


4 7'^ 


4 59 


4 7S 


4 59 


4 7S 


4 595 


0.95 


NLO 


4.77 


4.50 


4.78 


4.511 


4.79 


4.52 


0.55 


NNLO 


4.95 


4.54 


4.96 


4.55 


4.96 


4.55 


0.69 


/(soft = 2.5 GeV 

A^hard = 5 GcV 

A*fact = 5 GeV 


LO 


4.79 


4.63 


4.79 


4.63 


4.79 


4.63 


0.87 


NLO 


4.82 


4.505 


4.83 


4.517 


4.84 


4.525 


0.40 


NNLO 


5.08 


4.595 


5.08 


4.595 


5.09 


4.6 


0.69 


/(soft = 2 GcV 

yWhard = 5 GeV 

Aifact = 5 GeV 


LO 


4.84 


4.67 


4.84 


4.666 


4.84 


4.66 


0.72 


NLO 


4.84 


4.49 


4.85 


4.505 


4.86 


4.514 


0.24 


NNLO 


5.21 


4.66 


5.21 


4.66 


5.21 


4.66 


0.69 



TABLE L The kinetic mass my^^a for = 1 GeV as a function of the soft renormahzation 
scale for the moments n = 14, 16, 18 and successive approximations for Mia evaluated for NNLO 

mkin- 

used as an input. For comparison we quote here also the values of the M.ia moment for the 
pole mass, which corresponds to the formal limit /iq ^ in our approach. For //soft = 3.5 
GeV the NNLO pole mass equals to 4.95 GeV. Galculating M-u with this mass we obtain: 
M\2 = 0.2, M^l^ = 0.25 and M^f^^ = 0.69. 

Prom Table 1 one can see two things - the perturbation theory for the pole mass behaves 
in a way, that does not show any sign of convergence; the NNLO corrections normally exceed 
the NLO ones. Moreover, the pole mass strongly depends on the soft renormahzation scale. 
This picture is consistent with the expectation of the irreducible ambiguity of order Aqcd 
in the pole mass. 

For the low-scale mass the situation is different in both respects. The first terms of the 
perturbation series are sign alternating (if taken seriously, this feature signals that we are 
on the right way). Also, the low-scale running mass exhibits only a moderate dependence 
on /igoft in a relatively wide range of the soft renormahzation scale. The width of this range 
depends on the order of perturbation theory we consider; it also depends on the initial 
value of as{Mz). The most stable picture emerges at the NLO, while the inclusion of the 
NNLO effects makes the result less stable. This partial loss of stability is the consequence 
of the very large value of the second order correction 02 to the perturbative quark antiquark 
potential. 

If the sign alternating behavior of the perturbation series for mkin is taken seriously, 
one can try to perform some transformation of the perturbation series to accelerate the 
convergence. One of numerous options is the Euler transformation. It is interesting to 
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observe that this transformation indeed brings the values of mkin obtained for various values 
of A^soft closer to each other. 

The Euler transformation works in the following way. Imagine we have a series f{z) = 
l)"'c„z"'. Then the faster convergent approximation for f{z) is given by: 

/(.) ^ ^ (eo - (e. - c.) (^) + (c, - 2c, + c„) {^J + ...). (72) 

We restricted our consideration to two orders of perturbation theory which we can use to 
determine the expression for the mass. 

We then identify z in the previous formula with a(/isoft)- Then we use the results in 
Table 1 for n = 16 and fisoft = 4.5, 2.5, 2 GeV. For the sake of illustration, we present such 
calculation for yUsoft = 4.5 GeV: 



m = 4.57 - 0.07 + 0.02 = 4.57(1- 0.07a, + 0.09ay, a, = a,(4.5) = 0.22. (73) 
We then use Eq. (^2|) and the above formula to obtain a faster convergent series: 



4 57 / a o? \ 

mkin = — 1 + 0.93-^ + 0.95 . , ^ + ... = 4.49 GeV. (74) 

1 + as V 1 + as (1 + <^sY j 

The result does not change notably, as compared to the "naive" summation of the a, series, 
indicating that in this case the numerical value of a, is fairly small. 

In the same way, we obtain the new values of the kinetic masses for the yUsoft = 2.5 and 
2 GeV. The results are: for /igoft = 2.5 GeV, we obtain mkm = 4.52 GeV and for /^soft = 2 
GeV, mkin = 4.53 GeV. 

Clearly, such transformation cannot be rigorously justified; however, the fact that the 
numbers come closer to each other looks gratifying. Moreover, both numbers appear to 
become closer to the value of the NNLO mass, which can be obtained by examining the 
region of relative stability /isoft > 2.5 — 3 GeV, which we use below to estimate the value of 
the h quark mass. 

Let us also comment on how the choice of the value of a^iMz) is reflected on our result. 
The important point is that this dependence is rather moderate, since the change in the 
initial value of as{Mz) is roughly equivalent to the change in /igoft- As we always work in 
the region where the dependence of mkin on the normalization scale of the coupling constant 
is relatively weak, the same equally applies to the dependence on the initial value of a, at 
the Z-resonance. The dependence on the factorization scale is much weaker than on the 
soft renormalization scale and amounts to a variation of at most 20 MeV in value of the 
kinetic mass. 

We arrive finally at the following estimate for the kinetic mass extracted from the QCD 
sum rules: 

mkin(lGeV) = 4.56 ± 0.06 GeV. (75) 

We stress, that the error in the above estimate is primarily of the theoretical origin. The 
experimental errors in the masses and electronic decay widths of the T resonances, as well 
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as a poor knowledge of the continuum part of the observable spectrum, are minor effects as 
compared with, e.g., the soft scale dependence of the kinetic mass. For the above estimate, 
we use our results for /igoft between 4.5 and 2.5 GeV. For lower scales, the perturbation series 
for the moments do not look reliable enough. 

After the value of the pole mass is found, one can obtain the estimate of the MS mass 
rfi. To order 0{a1) the corresponding equation was given in UTG). One obtains: 



m(m) = 4.20 ±0.1 GeV. (76) 

Using Eq. (^, one can derive the evolution equation with respect to /ig and calculate 
the kinetic mass at different normalization scales. The only thing to be remembered is that 
a choice of the normahzation scale /iq is limited by two inequalities: Aqcd -C /xq <^ m. 
From this point of view, the pole mass, which formally corresponds to the limit /xg — > 0, 
is seen to be a completely artificial notion, since it includes a nonperturbative contribution 
treated in terms of the perturbation theory. 

Our NNLO results may look incomplete since we use relation (|^ between the pole and 
kinetic masses, which is valid to ^qoI order, while formally one needs to know this relation 
more accurately, to /xga^ order. However, a careful examination based on the BLM-estimate 
(D of the 0{iJiQa\) terms, shows that unknown corrections in that order in the relation 
between the pole mass and the kinetic mass cannot drastically change our results for the 
kinetic mass in the NNLO approximation. 

In any case, working with the values of the /^soft where perturbation theory seems to be 
reliable, we think we can provide a reasonable estimate for the value of the low-scale running 
mass to NNLO. 



VIII. CONCLUSIONS 



In this paper, we have determined the h quark low-scale running mass [|T^ from the 
analysis of the QCD sum rules in the next-to-next-to-leading order. We have shown that 
the use of this mass significantly improves the convergence of the perturbation series for the 
moments of the spectral density. As the result of our analysis we obtain the value of the 
kinetic mass normalized at 1 GeV: mkin(lGeV) = 4.56 ± 0.06 GeV and the corresponding 
value of the MS mass fn: fn{fn) = 4.20 ±0.1 GeV. 

In our opinion, the pole mass of the b quark cannot be reliably determined from the 
sum rule analysis. We have shown, that the NNLO order corrections to the pole mass are 
typically larger than the NLO ones. Also the value of the pole mass of the b quark is very 
sensitive to the scale of the strong coupling constant, that is used in the analysis. We think 
that these features are in accord with the fact that the pole mass of the quark cannot be 
reliably defined theoretically and suffers from an irreducible ambiguity of the order of Aqcd 

iiH. 



As a byproduct of our study we have obtained the NNLO analytical expression for the 
cross section e^e~ — QQ of the quark antiquark pair production in the threshold region. 
We have also given the NNLO expressions for the energy levels and the wave functions at 
the origin for the ^5*3 bound states of QQ. 
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